Flux growth and physical properties of Mo 3 Sb 7 single crystals 
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Millimeter sized single crystals of M03SD7 are grown using the self-flux technique and a thor- 
ough characterization of their structural, magnetic, thermal and transport properties is reported. 
The structure parameters for the high-temperature cubic phase and the low-temperature tetragonal 
phase were, for the first time, determined with neutron single crystal diffraction. Both X-ray pow- 
der diffraction and neutron single crystal diffraction at room temperature confirmed that M03SD7 
crystallizes in Ir3Ge7-type cubic structure with space group ImSm. The cubic-tetragonal structure 
transition at 53 K is verified by the peak splitting of (4 0) reflection observed by X-ray single crys- 
tal diffraction and the dramatic intensity change of (12 0) peak observed by neutron single crystal 
diffraction. The structural transition is accompanied by a sharp drop in magnetic susceptibility, 
electrical resistivity, and thermopower while cooling. A weak lambda anomaly was also observed 
around 53 K in the temperature dependence of specific heat and the entropy change across the tran- 
sition is estimated to be 1.80 J/molMo K. The temperature dependence of magnetic susceptibility 
was measured up to 750 K and it follows a Curie- Weiss behavior above room temperature. Analysis 
of the low-temperature magnetic susceptibility suggests a spin gap of 110 K around 53 K. A typi- 
cal phonon thermal conductivity was observed in the low temperature tetragonal phase. A glassy 
phonon thermal conductivity above 53 K suggests a structural instability in a wide temperature 
range. Superconductivity was observed at 2.35 K in the as-grown crystals and the dimensionless 
specific heat jump AC(T)/7 n T c was determined to be 1.49, which is slightly larger than the BCS 
value of 1.43 for the weak-coupling limit. 



PACS numbers: 74.25.-q, 74.25.fc,75.40.Cx 
INTRODUCTION 



Despite a relatively low superconducting transition 
temperature T c = 2.08 K,p] the Zintl compound M03SD7 
has attracted considerable interest due to the possible in- 
volvement of magnetism in superconducting pairing, [2]-[5] 
and promising thermoelectric performance with proper 
doping. [5H5] MoaSby is the only known compound in the 
Mo-Sb binary system. It crystallizes in a Ir3Ge7-type 
cubic structure with space group Im3m at room temper- 
ature. The Mo-sublattice is characterized by a three di- 
mension network of Mo-Mo dumbbells formed by nearest 
neighbours (NN) and octahedral cages at the body center 
positions formed by next nearest neighbours (NNN). The 
competition between strong antiferromagnetic NN and 
NNN interactions is proposed to drive a structure transi- 
tion at T* = 53 K from cubic to tetragonal (space group 
IA/mmm) to relieve spin frustration. |10ffT2"] This struc- 
ture transition further shortens the NN Mo-Mo distance 
along the crystallographic c-axis and forms spin singlet 
dimers. ^iSR, inelastic neutron scattering, and magnetic 
susceptibility measurements on polycrystalline samples 
suggest that the formation of a 120 K spin gap accom- 
panies the structural transition upon cooling. |2 H] With 
a structure transition and a spin gap above T c , Mo3Sb7 
is a unique system for exploring the interplay between 
magnetism, structure instability, and superconductivity. 



Despite much effort, the pairing mechanism at low 
temperatures and the correlation between the structural 
transition, magnetic frustration, and superconductivity 
are still controversial. Thorough studies of the intrinsic 
properties of Mo3Sb 7 and related compounds have been 
hampered by the unavailability of sizable single crystals. 
The first attempt to grow MosSb7 single crystals was 
performed by chemical transport reactions and resulted 
in small crystals adequate for determining the crystal 
structure. [H] Millimeter sized single crystals suitable for 
some physical property measurements were first grown 
by melting Sb in a thick Mo tube. [I] Upon cooling from 
1000°C to 700°C at the rate of 3°C/h, single crystals 
grow on the walls of the Mo tube. Crystals grown in this 
manner become superconducting at T c = 2.08 K, and en- 
abled several studies of the intrinsic properties. [T| HA\ [l"5] 
However, these crystals are still small in size. More im- 
portantly, it's difficult to grow crystals doped with other 
elements, especially substitutions at the Mo site. 

In this paper, we report the flux growth of millimeter 
sized Mo3Sb 7 single crystals and a thorough characteri- 
zation of the crystals by x-ray powder diffraction, x-ray 
and neutron single crystal diffraction, magnetic suscep- 
tibility, electrical resistivity, thermoelectric power, ther- 
mal conductivity, and specific heat measurements. Crys- 
tals grown by the self-flux technique posses the highest 
T c reported to date and are free of magnetic impuri- 
ties. Neutron single crystal diffraction was performed, for 
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the first time, to determine the structure parameters of 
the cubic phase at room temperature and the tetragonal 
phase at 4K. The structural transition is accompanied 
with anomalies in the temperature dependence of electri- 
cal resistivity, magnetization, specific heat, thermopower 
and thermal conductivity. 



EXPERIMENTAL DETAILS 

Single crystals of M03SD7 were grown out of Sb flux. 
The starting materials are elemental Mo powder (Alfa, 
99.999%) and Sb (Alfa, 99.9999%). Mo powder was first 
reduced in flowing Ar balanced with 4% H 2 for 12 hours 
at 1000°C. Reduced Mo powder and commercial Sb shot 
were mixed in the ratio of Mo:Sb = 1:49 and placed in 
a 2 ml AI2O3 crucible. A catch crucible of the same size 
containing quartz wool was mounted on top of the growth 
crucible and both were sealed in a silica ampoule under 
approximately 1/3 atmosphere of argon gas. The sealed 
ampoule was heated to 1000°C in 5 hours and homo- 
geneized for 12 hours before cooling to 700° C over 100 
hours in a programmable box furnace. At 700°C, the Sb 
flux was decanted from the MosSby crystals. 

Elemental analysis of the crystals was performed us- 
ing a Hitachi TM-3000 tabletop electron microscope 
equipped with a Bruker Quantax 70 energy dispersive 
x-ray system. X-ray diffraction from oriented single 
crystals and from powders ground from crystals was 
performed on PANalytical X'Pert Pro MPD powder 
x-ray diffractometer using Cu K Q i. The room tem- 
perature x-ray powder diffraction pattern was refined 
by the Rietveld method using Fullprof. [16] Single crys- 
tal neutron diffraction was measured at HB-3A four- 
circle diffractometer at the High Flux Isotope Reac- 
tor at the Oak Ridge National Laboratory. A neu- 
tron wavelength of 1.003 A was used with a bent per- 
fect Si-331 monochromator. [T?] The data were refined by 
the Rietveld method using Fullprof. Magnetic proper- 
ties were measured with a Quantum Design (QD) Mag- 
netic Properties Measurement System (MPMS) in the 
temperature interval 1.8 K < T < 750 K. The temperature 
dependent specific heat and electrical transport data 
were collected using a 14 Tesla QD Physical Proper- 
ties Measurement System (PPMS) in the temperature 
range of 1.9K<T<300K. Thermal conductivity and 
thermopower were measured from 1.9 K to 300 K using 
the Thermal Transport Option (TTO) from QD and a 
9 Tesla PPMS. A rectagular bar with the dimension of 
0.6 mm x 0.7mm x 6 mm was cut from a large crystal and 
used for the TTO measurement. Silver epoxy (H20E 
Epo-Tek) was utilized to provide mechanical and thermal 
contacts during the thermal transport measurements. 
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FIG. 1: (color online) Powder x-ray diffraction pattern of pul- 
verized M03SD7 single crystals recorded at room temperature. 
The observed (open circle) and calculated (solid line) profiles 
are shown on the top. The vertical marks in the middle are 
calculated positions of the Bragg peaks. The dashed line in 
bottom of the plot is the difference between calculated and 
observed intensities. Inset shows a piece of Mo3Sb7 single 
crystal. 

RESULTS AND DISCUSSIONS 
Crystal growth and phase stability 

As-grown crystals are normally cube-like with trun- 
cated corners as shown in the inset of Fig.[TJ The typ- 
ical dimension of a crystal is about 2-3 mm. Elemental 
analysis confirmed that the crystals are stoichiometric 
Mo3Sb7 with no observable deviation. Figure 1 shows the 
x-ray powder diffraction (XRD) data collected on ground 
Mo3Sb7 single crystals. The XRD data were analyzed by 
Rietveld refinement using the FullProf software and the 
Rietveld fit profiles are also presented in Fig.[T] The XRD 
data indicate that the crushed crystals are single phase 
Mo3Sb7. Weak reflections from Sb were observed in some 
measurements. This small amount of Sb is from residual 
flux on the surface of the crystals after the decanting and 
can be removed mechanically or by heating the crystals 
in evacuated and sealed quartz tube at 600°C with the 
cold end at room temperature. 

Growths with various charge/flux ratios were per- 
formed in order to optimize the growth and the yield. 
We noticed that raising the charge to flux ratio or cooling 
faster leads to smaller Mo3Sb7 crystals. No sizeable crys- 
tals were obtained using the above process once the Mo 
content is over 5 at%, even with a 24-hour homogeniza- 
tion at 1150°C before cooling. From the presently avail- 
able phase diagram. [TH] Mo3Sb7 decomposes at 780°C 
into (Mo) and liquid (Sb). However, a test growth 
stopped at 800°C also yields Mo3Sb7 crystals from a 
Mo:Sb=l:49 melt. This suggests that there is a wide 
temperature and composition range where Mo3Sb7 and 
liquid (Sb) coexist. However, the liquidus temperature 
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FIG. 2: (color online) Plot of the x-ray diffraction pattern for 
(400) cll bic reflection in Mo3Sb7 as a function of temperature 
measured on a piece of oriented single crystal. The stronger 
intensity of (004)tet than that of (400)tet is due to the domain 
distribution. 



should increase rapidly with increasing Mo content in 
the melt, which limits the charge/flux ratio during crys- 
tal growth. A Mo/Sb ratio of 1:49 has been found to 
reproducibly grow the largest crystals. 



Structure 

The Rietveld refinement of the x-ray powder data 
shown in Fig.[l] confirmed that as-grown MoaSby single 
crystals crystallize in the Ir3Ge7 type structure with the 
cubic Im3m (space group 229) symmetry at room tem- 
perature. The room temperature structure is further con- 
firmed by a single crystal neutron diffraction measure- 
ment. Table 1 shows the crystallographic parameters 
refined for the room temperature x-ray powder diffrac- 
tion pattern and the neutron single crystal diffraction 
data. The results agree well with each other and previ- 
ous reports, [mn^] 

A structural transition from cubic to tetragonal was 
reported to take place at 53 K in polycrystalline samples. 
Single crystal x-ray and neutron diffraction were used to 
study this transition in our crystal. Figure [2] shows the 
evolution of the cubic (4 0) peak with temperature from 
11 K to 100 K. Below 54 K, the peak splits signaling the 
lowering of the symmetry. At 11 K, the a/c ratio was 
obtained to be 1.002, consistent with previous reports 
from powder x-ray diffraction. 1 1 2] 

Figure [3] shows the temperature dependence of the in- 
tegrated intensity of the cubic (12 0) nuclear peak mea- 
sured with neutron single crystal diffraction. Two fea- 
tures are noteworthy: (1) the (12 0) broadens below 
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FIG. 3: (color online) Temperature dependence of the inte- 
grated intensity of the cubic (12 0) reflection measured by 
neutron single crystal diffraction on a Mo3Sb7 single crystal. 
Inset shows the peak splitting and the fitting of low temper- 
ature peaks. 



T* = 53 K signaling the structure transition. Using the 
peak width of (12 0) fitted by Gaussian function at high 
temperature, the peak splitting was obtained by fitting 
with two Gaussian functions (see inset of Fig. [3]). (2) a 
rapid increase of the (12 0) peak intensity was observed 
upon cooling through the structural transition. At such 
a high q, the magnetic contribution to the intensity, if 
any, should be negligible. This peak intensity change is 
due to the weakened extinction associated with the cubic- 
tetragonal structural transition at low temperatures. Ta- 
ble 2 shows the crystallographic parameters refined for 
the neutron single crystal diffraction data collected at 
4 K. Refinement was performed with the space group 
I4/mmm (space group 139) and the crystal twinning was 
considered in the refinement. We noticed that the atomic 
displacement parameters (ADPs) at room temperature is 
3-5 times larger than those at 4K. This change of ADPs 
also contributes to the intensity change of (12 0) peak; 
however, a large decrease in extinction in this crystal 
occurs at the structure transition. The larger ADPs at 
300 K may not be induced by the structure instabilities 
in the cubic phase. In RU3S117 with the same structure 
but no structural instabilities, similar ADPs were also 
observed. [TS] 

Figure [4] shows the crystallographic structure with the 
parameters listed in Table 1 and Table 2. In the cubic 
phase above 53 K, the Mo sublattice could be looked as 
three orthogonal -Mo-4Mo-Mo- chains as shown in Fig|4] 
(a) and (c). In the tetragonal phase, Mo-Mo nearest 
neighbors along c-axis dimerize. As illustrated in Figj4] 
(d) and (e), the dimerization distinguishes the -Mo-4Mo- 
Mo- chains along c-axis from those along a- and b-axes. 
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TABLE I: Structure parameters for M03SD7 with the space group Im3m (no. 229) at room temperature. Atomic coordinates 
are : Mo (x,0,0), Sbl (0.25,0,0.5), and Sb2 (x,x,x). X-ray data were obtained from Rietveld refinement of powder diffraction 
pattern collected on pulverized single crystals. Neutron data were measured on a piece of single crystal of 30 mg at HB-3A. 
Lattice parameter a was determined to be 9.569(2)A by x-ray and 9.543(2)A by neutron diffraction. 



atom 


site 


z(Xray/Neutron:X/N) 


t/ lso (A 2 )(X/N) 


Mo 


12e 


0.3438(3)/0.34340(16) 


0.0060(3)/0.0041(5) 


Sbl 


12d 


1/4 


0.0084(3)/0.0068(5) 


Sb2 


16/ 


0.16268(9)/0.16219(10) 


0.0081(3)/0.0055(4) 



TABLE II: Structure parameters for Mo3Sb7 with the space group IA/mmm (no. 139) at 4 K . Neutron data were measured 
on a piece of single crystal of 30 mg at HB-3A. The lattice parameters are a = 9.551(4)A, and c = 9.523(2)A with a,/c — 1.003. 
Rf=0.0373, x 2 = 3.563 for 210 F >2(j(F ). 



atom site x y z Ui SO (A 2 ) 

MoT 81 0.3436(2) 0.0016(6) 

Mo2 4e 0.3442(4) 0.0012(9) 

Sbl Ad 1/2 1/4 0.0013(11) 

Sb2 87 0.2501(6) 1/2 0.0012(6) 

Sb3 16m 0.16232(13) 0.16232(13) 0.1621(3) 0.0010(4) 



Mo2 




FIG. 4: (color online) Crystal structure of MosSb7 at (a) 
room temperature, and (b) 4K. The -Mo-4Mo-Mo- structure 
unit and selected bond lengths (in unit of A) at (c) 300 K, (d) 
4K along a-axis, and (e) 4K along c-axis. 



Superconducting transition 

Superconductivity is observed to occur below 
T c = 2.35 K, as indicated by the magnetization, electrical 
resistivity, and specific heat data shown in Fig. [5] 
Magnetization was measured in both zero-field-cooling 
(ZFC) and field-cooling (FC) modes with an applied 



magnetic field of 10 Oc. The observed diamagnetism 
shown in Fig. [5ja) confirms bulk superconductivity. The 
onset of the superconducting transition temperature is 
^2.35 K, which agrees well with that determined from 
the temperature dependence of resistivity (also shown 
in Fig[5]ja)). Figure [5jb) shows the low temperature 
specific heat data plotted as C p /T vs. T 2 . A specific 
heat jump at 2.35 ± 0.05 K was well resolved at the 
onset of superconductivity. High magnetic fields (>2T) 
suppress superconductivity and the specific heat jump 
disappears as shown in Fig.[5^b). Literature values of 
T c vary from 2.02 K to 2.30 K, with the variation likely 
coming from sample quality. [2 [TU HS1 I2H] The current 
as-grown crystals posses the highest T c reported to date. 

The fitting of specific heat data in the range 
8 < T 2 < 60 K 2 to the standard power law, C p = 7T + /3T 3 
yields 7 = 31.8(2) mJ/molK 2 and /3 = 0.64(1) mJ/molK 4 , 
where 7 is the Sommerfeld electronic specific heat co- 
efficient and ft the coefficient of the Debye T 3 lat- 
tice heat capacity at low temperatures. The latter 
gives the Debye temperature 9d with the following re- 
lation 0jj = (127r 4 N J 4ksn/5/?) 1 / 3 , where n is the num- 
ber of atoms per formula unit, is Avogadro's con- 
stant and ke is Boltzmann's constant. With n = 10 
and j3 = 0.64(1) mJ/molK 4 , the Debye temperature is 
#d = 248K for Mo 3 Sb 7 . We note that the Debye tem- 
perature is similar to that reported by Candolfi et al.,[2] 
but is smaller than that reported by Tran et al. [3], [21] 
The dimensionless specific heat jump AC(T)/7„T C was 
determined to be 1.49, close to 1.56 reported by Tran et 
al.[21 This value is slightly larger than the BCS value of 
1.43 for the weak-coupling limit. 

The density of states at the Fermi level D(Ej?) can be 
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FIG. 5: (color online) (a) The temperature dependence of 
magnetic susceptibility (left) measured in both zero-field- 
cooling (ZFC) and field-cooling (FC) modes with a field 
of 10 Oe and electrical resistivity (right). (b) The low- 
temperature specific heat of Mo3Sb7 single crystal divided 
by temperature as a function of temperature square mea- 
sured under T (open circles) and 12 T (solid circles) mag- 
netic fields. The solid lines are guide to the eye. The dashed 
line is the linear fitting of the low temperature specific heat 
as described in the text. 



obtained from the Sommerfeld coefficient 7 according to: 



7 : 



7r 2 fc| 



D(E F )(l + X e - ph ) 



(1) 



where A e _ p h is the electron-phonon coupling constant 
and could be estimated with the McMillan equation: 



A 



e—ph 



1.04 + ij,*ln(6 D /lA5T c ) 
(1 - O.62/x*)Zra(0 D /1.45r c ) - 1.04 



(2) 



where zx* is the Coulomb pseudopotential and usually 
taken between 0.1 and 0.15. With fi* =0.15, 9d = 248K, 
and T c = 2.35 K, X e - P h is 0.59 consistent with a previous 
report. [21] The density of states at the Fermi level D(Ep) 
is calculated to be 8.5states/eVf.u., where f.u. means 
formula unit. The value is comparable to that obtained 
by theoretical calculations. [22 




FIG. 6: (color online) The temperature dependence of mag- 
netic susceptibility of Mo3Sb7 single crystal measured un- 
der a magnetic field of 60kOe in the temperature range 
2 K < T < 750 K. The solid curve below 50 K shows the fitting 
as described in text. The lower inset shows the temperature 
derivative of magnetic susceptibility, d^(T)/dT. The upper 
inset shows the linear fitting of x(T) _1 . 



Physical properties in normal state 

Figure[6] shows the temperature dependence of the 
magnetic susceptibility, x(T), measured with an applied 
field of 60 kOe in the temperature range 2 K < T < 750 K. 
For x{T) measured above room temperature, the con- 
tribution from the sample holder was appropriately cor- 
rected. With decreasing temperature, x(T) increases and 
shows a broad maximum around T m ax — 175 K which is 
followed by a sharp drop around 50 K. The sharp drop is 
best illustrated by the derivative dx/dT, which shows a 
distinct maximum around 50 K as shown in inset of Fig. 5. 
The features below room temperature agree with previ- 
ous reports. [H [3] It's noteworthy that below 20 K x(T) 
shows negligible temperature dependence. While in all 
previous reports on the magnetic properties of polycrys- 
talline samples, a Curie- Weiss (CW)Tike tail is observed 
below 50 K, which has been attributed to unpaired spins. 
The absence of the CW-like tail further confirms the high 
quality of our crystals. The isothermal magnetization 
(not shown) measured at 5 K and room temperature in- 
dicates no ferromagnetic component. 

As shown in Fig.[6j x(T) below 50 K was fit with a 
spin gap function x(T) = x(0) + &exp(-A/keT), where 
x(0) is a temperature independent term, b is a con- 
stant, A is the spin gap, ks is the Boltzmann constant. 
The best fitting was obtained with the following param- 
eters: x(0) = 3.7xl0~ 4 cm 3 /mol, b = 5.8 x 10" 4 cm 3 /mol, 
and A/ks = 110(6) K. The gap is similar to that reported 
on polycrystalline samples. [3] 

The magnetic susceptibility of Mo3Sb7 in the temper- 
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ature interval 220 K < T < 300 K has been reported [3] to 
follow Curie- Weiss law. The upper inset of Fig. [6] shows 
the linear fitting of x(^)" 1 above 300 K. The fitting yields 
# = -1078K and p e ff = 1.67 /is /Mo. The effective mag- 
netic moment is close to the expected value for Mo ions 
with S=l/2. While such a large negative Weiss constant 
suggests strong antiferromagnetic interaction between lo- 
calized Mo moments. 

Figure [7] shows the temperature dependence of elec- 
trical resistivity in the temperature range 3 K ~ 300 K. 
The electrical resistivity at room temperature is about 
160/if2.cm. As highlighted by the solid line in the Fig- 
ure, at temperatures above ^150 K, the electrical resis- 
tivity shows a linear temperature dependence consistent 
with the report by Bukowski et al.pQ Starting around 
50 K, the electrical resistivity shows a sharp drop. This 
sharp drop can be better illustrated by the derivative 
of the electrical resistivity shown in the upper inset. 
At even lower temperatures, the resistivity levels off 
with cooling and finally drops to zero at T c . The low 
temperature electrical resistivity in polycrystalline sam- 
ples has been found to fit a p(T) = po + &T™ law with 
n = 2. 2] The dashed curve in the lower inset shows the 
fitting, and obviously the fit is not acceptable for the 
data below ~40K. A much better fit can be obtained 
with n around 2.4, which has been observed when fitting 
p(T) measured on crystals from different batches. We 
also tried to fit the low temperature electrical resistivity 
with a gap function p(T) — p + cT + dexp(-Ap/KsT) 
which has been used to argue against the spin fluctu- 
ation scenario and point to a gap feature. |3 As shown 
in the lower inset of Fig.[7J a good fit could be ob- 
tained with po = 91.5 /i.fi.cm, c = 5.58 XlO^fi.cm.K -1 , 
d = 9.36 x lO^flcm, and Ap/K B = 107(2) K. 

The superconducting pairing mechanism at low tem- 
peratures is still under hot debate despite much effort. 
Nazir et al. compared the electronic properties of the 
cubic and tetragonal phase in MosSby by density func- 
tional theory calculations. [35] A higher density of state 
at Fermi level is found in the tetragonal phase which 
has a higher propensity to magnetism than the cubic 
phase. Spin fluctuations were argued to coexist with su- 
perconductivity in Mo3Sb7.[2] One important evidence 
for this argument is the quadratic temperature depen- 
dence of electrical resistivity and magnetic susceptibility 
at low temperatures in the tetragonal phase. The ab- 
sence of the expected quadratic temperature dependence 
signals that the electron-paramagnon interaction should 
be moderate. [5] 

X-ray diffraction on oriented single crystals (Fig. 2) 
and neutron single crystal diffraction (Fig 3) clearly 
show the symmetry lowering from cubic to tetragonal 
while cooling. The shortened Mo-Mo bond length along 
the crystallographic c-axis signals only 1/3 Mo-Mo pairs 
dimerize. The fact that spin singlet state is only a frac- 
tion of the magnetic ground state explains the large mag- 




UK) 



FIG. 7: (color online) The temperature dependence of elec- 
trical resistivity of Mo3Sb7 single crystal. The solid line high- 
lights the linear temperature dependence of electrical resistiv- 
ity at high temperatures. The upper inset shows the deriva- 
tive of electrical resistivity, dp/dT. The lower inset compares 
the fitting with power law and a gap function as described in 
text. 



netization at low temperatures. On the other hand, 
the above feature distinguishes Mo 3 Sb7 from other low 
temperature superconductors by the coexistence of spin 
dimer singlet and superconductivivty. In a spin dimcr 
system, spin fluctuations from the triplet excitations 
could induce an attractive interaction between conduc- 
tion electrons leading to the superconducting pairing. [23] 
An preliminary study on Te or Ru doped Mo3Sb 7 sug- 
gests that a few percent of Te or Ru substitution elimi- 
nates the structure transition and superconductivity sur- 
vives in the cubic phase. Thus the structure transition 
or possible spin fluctuations from the triplet excitation 
cannot be the only driving force for the superconducting 
pairing. 

Figure [8] shows the temperature dependence of the spe- 
cific heat, Cp(T), measured under 0T and 12 T applied 
magnetic fields. The magnetic field of 12 T suppresses 
superconductivity and the specific heat jump at T c dis- 
appears and the details were shown in Fig.[5|h). The 
heat capacity at 200 K attains a value of 226J/mol K, 
which is a little smaller than the classical high tempera- 
ture Dulong- Petit value of 3nR = 249 J/mol K at constant 
volume, where R is the molar gas constant and n = 10 
is the number of atoms per formula unit. This agrees 
with the determined Debye temperature, which is slightly 
below room temperature. The inset of Fig. [8] highlights 
the details around 50 K. A weak lambda anomaly is well 
resolved around T = 53K. No field dependence was ob- 
served for this anomaly with magnetic fields up to 12 T. 

In order to estimate the entropy change around 53 K, 
we subtracted the lattice contribution to the specific heat 
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FIG. 8: (color online) (a) The temperature dependence of 
specific heat of M03SD7 single crystal measured under and 
12 T magnetic fields. The inset highlights the weak lambda 
anomaly around 53 K. (b) C p /T versus T and the numeri- 
cal integration of the entropy S after subtracting the lattice 
specific heat using Ru 3 Sn 7 as phonon reference. 



using Ru 3 Sn 7 as a reference. [2~4] Figure ^b) shows the 
residual specific heat plotted as AC p /T vs T. The en- 
tropy change, which is calculated as J AC P /T dT, is 
1.80 J/molMoK and about 1/3 of the expected value. 
This suggests that the entropy change is from Mo spins. 
As illustrated in Fig.[8]jb), about 60% entropy change 
takes place above 53 K. 

Figure [9] shows the evolution of the thermopower, 
a(T), with temperature. A room temperature value of 
~16//V/K agrees with the measurement on a polycrys- 
talline sample. [2"D] Above 100 K, cv(T) shows a nearly lin- 
ear temperature dependence. Below about 60 K, a(T) 
decreases quickly when cooling; this sharp drop resem- 
bles that in the temperature dependence of electrical 
resistivity around 53 K. Below 30 K, ev(T) stays around 
zero. It's interesting to note that p(T) levels off and is 
nearly temperature independent below 30 K. The ther- 
mopower is proportional to the logarithmic derivative of 
the density of states (DOS) with respective to energy at 
the Fermi level, while the electrical conductivity increases 
with DOS. The observed change of a(T) and p(T) sug- 
gests an electronic structure change across the structure 
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FIG. 9: (color online) The temperature dependence of ther- 
mopower of Mo3Sb7 single crystal. The dashed line is zero 
line. 

transition. 

Figure [10] shows the temperature dependence of the 
thermal conductivity in the normal state. With decreas- 
ing temperature, the thermal conductivity decreases from 
room temperature down to 50 K; below 50 K, a maximum 
near 20 K was observed. Since no long range magnetic 
order was observed in Mo 3 Sb7, possible heat carriers are 
electrons and phonons. Thus the total thermal conduc- 
tivity can be described as: K to tai = K e + K p , where K e is 
the electronic thermal conductivity, and k p the lattice 
thermal conductivity. n e was estimated from the elec- 
trical resistivity data using the Wiedemann-Franz law: 
K e =LT/p, where L is the Lorenz constant taken to be 
equal to 2.44 x 10~ 8 V 2 /K 2 , T is the absolute tempera- 
ture, and p is the electrical resistivity. Figure 9 shows 
the temperature dependence of K e and k p ; the latter was 
obtained by subtracting K e from the total thermal con- 
ductivity. K e decreases while cooling in the whole tem- 
perature range studied. We noticed that k p is rather low 
in the whole temperature range. This might be due to 
large number of heavy atoms in the primitive unit cell or 
strong electron-phonon scattering. The maximum near 
20 K for IS typical for a crystalline material though the 
maximum occurs at a low value. The striking feature for 
k p is a glassy behavior above 53 K where the structure 
transition takes place. This kind of phonon glass be- 
havior has been observed in a wide variety of materials, 
such as RV0 3 ,[25] RCo0 3 ,[2S] La 4 Ru 2 Oio, 27 and filled 
skutterudite antimonides.[5S] In the oxides with local- 
ized electrons, bond-length fluctuations induced by spin 
and/or orbital fluctuations disturb the phonon thermal 
conductivity through spin-orbital-lattice coupling. While 
in filled skutterudite antimonides, the rattling mode se- 
riously damps the lattice vibration. We note that neu- 
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FIG. 10: (color online) The temperature dependence of ther- 
mal conductivity of M03SD7 single crystal. See the text for 
the calculation of electron and phonon thermal conductivities. 

tron single crystal diffraction didn't observe abnormal 
atomic displacement parameters for either Mo or Sb in 
Moc.Sb7.f29] Thus, the phonon glass behavior of Mo3Sb7 
might signal the effect of the structural instability result- 
ing from the competing NN and NNN magnetic interac- 
tions. The cubic-tetragonal structural transition relieves 
the magnetic frustration thereby recovering the phonon 
heat transport in the tetragonal phase. This might be 
further enhanced by the reduced spin scattering due to 
the accompanied spin gap opening. 



SUMMARY 

Mo 3 Sb7 single crystals were successfully grown out of 
Sb flux. Our growth study suggests that there is a wide 
temperature and composition range where Mo3Sb7 and 
liquid (Sb) coexist. The sizeable crystals enabled the 
first neutron single crystal diffraction study of the high 
temperature cubic phase and the low temperature tetrag- 
onal phase. The physical properties are generally in line 
with previous reports on polycrystalline samples. How- 
ever, the following features are noteworthy: (1) Magne- 
tization, electrical resistivity and specific heat measure- 
ments observed superconductivity at 2.35 ± 0.05 K in the 
as-grown crystals, which is the highest T c reported to 
date for Mo3Sb7. The corresponding dimensionless spe- 
cific heat jump AC(T)/7„T C was determined to be 1.49, 
which is slightly larger than the BCS value of 1.43 for the 
weak-coupling limit. (2) A weak lambda anomaly was 
observed at 53 K in the temperature dependence of spe- 



cific heat and no field effect was observed with magnetic 
field up to 12 T. The entropy change across the transi- 
tion is estimated to be 1.80 J/molMoK and 60% entropy 
change takes place above 53 K. (3) The magnetic suscep- 
tibility shows negligible temperature dependence below 
^20 K suggesting the crystals are free of magnetic im- 
purities. Above room temperature, the magnetic suscep- 
tibility follows CW law and a moment of 1.67 /ig/Mo is 
close to the expected value for Mo ions with spin S = 1/2. 
Thus the electronic state of Mo has both localized and 
itinerant characters. (4) The electrical resistivity below 
50 K doesn't show a quadratic temperature dependence. 
Analysis of the transport and magnetization data sug- 
gests the opening of a spin gap of 110 K accompanying 
the structure transition. (5) A phonon glass behavior 
was observed above 53 K suggesting strong scattering 
of phonons by lattice instabilities. The abrupt change 
around 53 K in the temperature dependence of electrical 
resistivity and thermopower suggests an electronic struc- 
ture change across the structural transition. 

The sizeable high quality single crystals allow stud- 
ies of intrinsic properties of Mo 3 Sb7. Moreover, prelimi- 
nary growths have demonstrated that crystals with var- 
ious substitutional dopants, for either Mo or Sb, can be 
grown with controlled doping using the reported growth 
protocol. This enables various future studies that will ex- 
plore the relationship between the structure instability, 
magnetism and superconductivity, and the origin of the 
promising thermoelectric performance of doped Mo3Sb7. 
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